SOME ORTHOGONAL VERY- WELL-POISED s'/'t-FUNCTIONS 
THAT GENERALIZE ASKEY-WILSON POLYNOMIALS 



SERGEI K. SUSLOV 



Abstract. In a recent paper Ismail, Masson, and Suslov |Q have established a continuous orthog- 
onality relation and some other properties of a 2'/3i-Bessel function on a q-quadratic grid. Dick Askey 
suggested that the "Bessel-type orthogonality" found in at the 2¥'i-level has really a general 
character and can be extended up to the 8(^7-level. Vcry-well-poised g^Jy-functions are known as a 
nonterminating version of the classical Askey- Wilson polynomials ||^. Askey 's congecture has 
been proved in [Q. In the present paper which is an extended version of |33| we discuss in details 
some properties of the orthogonal 8</57-functions. Another type of the orthogonality relation for a 
very- well-poised sV'T-function was recently found by Askey, Rahman, and Suslov H]. 



1. Introduction 



The Askey- Wilson polynomials ^ are 

Pn{x) = pn{x;a,b,c,d) 

a"" {ab, ac, ad; g)„ 4V93 | ; q, q 

ab, ac, ad 



(1.1) 



where x = cos 6. These polynomials are the most general known classical orthogonal polynomials 
i, 0, i, PI, @, and H). 



see 



The symbol 4^93 in (1.1) is a special case of basic hypergeometric series W2 



/ 



riPsit) := 



E 



ai, 02, • • • , Cf-r 



\ 61,62, • • • ,bs 

(«!, 02, . . . , ar] q)r 



q, t 



[1.2) 



"0 {q,Kh,--- ,6,;g)„ 
The standard notations for the g-shifted factorials are 

n 



((-1 



k=l 



(oi, 02, . . . , ar, q)n ■= n ("fc! ^)"' 



:i.3) 
:i.4) 



k=l 
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and 

(a; g)oo := lim (a; (1.5) 

n—i-oo 

r 

(ai, 02, . . . , ar] q)oo ■= Yi i^^k] q)oo (1-6) 

k=l 

provided |g| < 1. For an excellent account on the theory of basic hypergeometric series see [p!2| . 
We shall just mention here that the s+i(/9s-series is called balanced if qaia2 ■ ■ ■ cis+i = &1&2 ■ ■ - bs and 
t = g; it is a very-well-poised if 0261 = 0362 = • . . = cts+i&s = qai, and 02 = q^Ja{-, 03 = —q^Jal- 

Askey and Wilson found the orthogonality relation 

f-TT p„(cos 6*; a, 6, c, d) Pm{cos 6; a, 6, c, (^e^*^, e"^*^; 



5, 



(g, a6, ac, ad, be, bd, cd; q)oo 
277 {abed; q)oc 
(1 — abcdq~^) {q, ab, ac, ad, be, bd, cd; q)^ 



de (1.7) 



{1 — abcdq^"^ ^) {abcdq ^;g)„ 
In the fundamental memoir |0, they studied in details many other properties of these polynomials. 

As is well known, the Askey-Wilson polynomials and their special and limiting cases are the 
simplest and the most important orthogonal solutions of a difference equation of hypergeometric 
type on nonuniform lattices (see, for example, 0, [0, §], 0, |jT2|, |2^, and [0). Recently Ismail, 
Masson, and Suslov have found another type of orthogonal solutions of this difference equation 
||l^-|T3. They considered the 2V^i-function, 

J,{z,r) = J,{x{z),r\q) (1.8) 
'ry (g^'+S -rV4;g), 



27 (g;g)c 



2^1 



as a g-analog on a g-quadratic grid of the Bessel function p 



Ismail, Masson, and Suslov established the following orthogonality property for the g-Bessel func- 
tion, 

J,{cose,r)J,{cose,r')- ^-- — — (1.10) 

^ ^ ^ ^ (g"e'^g'^e-*^gi-°e*^gi-°e-*^;g)^ ^ ^ 



' 00 



if r 7^ r'. 



(q^''+'^/^e'',q^''+'^/^e-'',q^''+'^/^e'',q^''+'^/^e-'';qY' dO = 



.2 (e2*^e-2^^;g)^ 



/oo 



V / 00 
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-47r(l - dJ^{x{a),r) VX(x(a),r) 



' oo 



if r = r'. Here r and r' are two roots of the equation 

Juia,r) = Jy{a,r') = 0, 



Vx(a) 



:i.l2) 



and J^{z,r) = Ji,{x{z),r), x{z) = - (q^ + q (x = cos 6', if = e*^), Re > —1, and < Re a < 

1. (See ||16|-||l^ for more details.) This is a g-version of the orthogonahty relation for the classical 
Bessel function 



X Jyirx) Jyir'x) dx 



if r 7^ r', 

]- {Jy+i{r)f if r = r', 



(1.13) 



under the conditions Jv{r) = Jv{r') = p^. 



Another example of an orthogonality relation of this type has been discovered recently. The 
following functions 



C{x) ■.=Cq{x; u) 



(1.14) 



[-uj^;q^)o^ 
-qu'^-^q^)^ 



2^1 



-qe^'', -qe-^'' 



2 2 

; q , -UJ 



and 



S{x) ■. = Sq{x] u) 

_ {-u;';q')^ 2q^l' 

{-quj^;q'^)^ 1-q 

X COS^ 2<^l 



:i.l5) 



2 2 

q , -u 



were discussed in [§] and as analogs of coscux and smujx on a g-quadratic lattice, respec- 
tively. Ismail and Zhang p2| expanded the corresponding basic exponential function in terms of 
"g-spherical harmonics" and later, together with Rahman [^, they extended this g-analog of the 
expansion formula of the plane wave from g-ultraspherical polynomials to continuous g-Jacobi poly- 
nomials. "Addition" theorems for the basic trigonometric functions (1.14)-(1.15) where found in 



3^ . Bustoz and Suslov |Tl| have established an orthogonality property. 



(:7(cos^;c<;)C(cos^;cj') 
^(cos 0] oj) S'(cos 6; ui') 
C{cos9] u) S{cos9; u') 



^-2^e. q 



(gi/2e2ie^gi/2e 


-2ie. 


( p2i9 -2ie. 


") 

/ oo 




(gi/2e2i9, gi/2e- 




( 2ie p-2ie 


/ oo 


1 c , c 



^qi/2^2ie^qi/2^-2ie.qy 



dd = 0, 



de = 0, 



de = 0, 



(1.16) 
(1.17) 
(1.18) 
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and 



vr 



Here uj and tu' are different solutions of the equation 



S 



Q(g^/^ + g-^/^);^)=o. (1.20) 

Bustoz and Suslov introduced the corresponding g-Fourier series and established several important 
facts about the basic trigonometric system and the g-Fourier series. 

Dick Askey has suggested that the orthogonality relation (l.lO)-(l.ll) can be extended to the 
level of very- well-poised gy^y-f unctions. The author was able to prove his conjecture in ||3^. Our 
main objective in this paper is to study the new orthogonality property for g^^y-function in details. 
This article is an extended version of the short paper |Q, originally submitted as a Letter, we 
include proofs of the main results etablished in to make this work as self-contained as possible. 



The paper is organized as follows. In Section 2 we consider difference equation of hypergeometric 
type on a g-quadratic grid and discuss a solution of this equation in terms of a very-well-poised 
gyjy-function. In the next two sections we derive a continuous orthogonality property for this sV7~ 
function. Section 5 is devoited to the investigation of zeros and asymptotics of this function and 
in Section 6 we evaluate the normalization constants in the orthogonality relation for the gyjy. In 
Section 7 we find an analog of the Wronskian determinant for solutions of difference equation of 
hypergeometric type. Some special and limiting cases of our orthogonality relation are discussed in 
Section 8. We close this paper by estimating the number of zeros on the basis of Jensen's theorem 
in Section 9. 



2. Difference Equation and Its gv^y-SoLUTioNS 
Let us consider a difference equation of hypergeometric type 




(2.1) 

on a g-quadratic lattice x{z) = - (q^ + g"^) with Xi{z) = x (^z + and Af{z) = Vf{z + 1) = 
f{z + l) — f{z). Here, in the most general case, 

a{z) = q-'^ (g^ - a) (g^ - b) (g^ - c) (g^ - d) , (2.2) 
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TiZ 



a[-z) - a[z) 
Vxi{z) 



1-q 



{abc + abd + acd + bed — a — b — c — d + 2(1 — abed) x) 



A = A. 



I — q (1 — abcdq 



Equation (2.1) can also be rewritten in self-adjoint form, 



A 



Vxi(^) 



(j{z)p{z^ 



Vu{z) 
S/xiz) 



+ \p{z) u{z) 



(2.3) 
(2.4) 

(2.5) 
(2.6) 



where p{z) is a solution of the Pearson equation, 

/\{a{z)p{z)) = T{z) Wxi{z). 

See p5| and |^ for details. 

Methods of solving of the difference equation of hypergeometric type (2.1) were discussed in [||, 
||13|| , and [^. As is well known, there are different kinds of solutions of this equation. For integer 
values of the parameter z/ = n = 0, 1, 2, . . . , the famous solutions of (2.1) are the Askey-Wilson 
4(y93-polynomials (see, for example, 0, ||^, and [jT^). For arbitrary values of this parameter 
solutions of (2.1) can be written in terms of 8 V'r-f unctions |18|, and Let us choose the 

following solution, Uy{z) = {x{z)] a, b, c; d), such that 



My(a;; a, b, c; d) 

_ {qa/d, beq", q^-^+'/d, q^~^~'/d] q)^ 
~ {q^-'a/d, be,q^+-/d,q^--/d; q)^ 
I 



(2.7) 



X sV^r 



aq 



d 



laq 



d 



d 



bd ed 



q, beq" 



laq 



\ 



d 



laq'" aq , q^-''^' q^-"" 
——, — -, ab, ae, — - — , — - — 
d d d d 



(beq^q'-'^/ad; q)^_ 



( q ^,abedq^ ^,aq^,aq 



+ 



{be,q/ad; q)^ 
{q~'^, abedq"'^, qb/d, qe/d; q)^ 



q, q 



(2. 



V 



ab, ae, ad 



{ab, ae, be, ad/q; q)^ 
{aq\aq-^;q)^ 



iq^+^/d,q^-'/d; g), 



4^3 



( q^-^/ad, beq", q^+'/d, q^-'/d 
\ qb/d,qe/d,q^ /ad 



q, q 



We have used Bailey's formula (III. 36) of |T2[ to transform (2.7) to (2.8). Both 4(^3-functions in 
(2.8) are balanced and converge when |g| < 1. The sV^r-function in (2.7) has a very-well-poised 
structure, it converges when \beq'^\ < 1; (2.8) provides an analytic continuation of (2.7). This s^^t 
can be transformed to another equivalent form by (III.23)-(III.24) of 
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A similar function was earlier discussed by Rahman ||2^ and in a recent paper he has found 
a g-extension of a product formula of Watson involving this function. Rahman has also shown that 

/ 1 — a;' 



2 

Thus, the sV^y-function in (2.7) can be thought as a g-extension of the hypergeometric function of 
Gauss. 

It is easy to see that Uu{x]a,b,c]d) defined by (2.7)-(2.8) is a function in x = {q^ + q~^) /2. 
Indeed, by (1.3) and (1.4), 

n-l 
" k=0 

where ^ = a, q/d, q^~'^/d and n = 1, 2, 3, . . . or oo. It is also important to mention that 



n-l 

n (l - {q'"" + abcdq"-^) q^ + ahcd g^^"^ 

fc=0 

n-l/ / _^ (l_g)2 



n ( 1 - ( 1 + ahcdq-^ - A J g*^ + ahcd g^^"^ 



fc=0 



4g3/2 



and 



{hc(f "/ad] 



n (1 - i fl + ^-i*- - ^ A + 5 9"*'^ 



ad \ ^q^f^ J ad 



by (1.3)-(1.4) and (2.4), where n = 1, 2, 3, . . . or 00. Therefore, Uy(a;; a, b, c; d) is, really, a function 
of At, as well. 



The definition of the u^, in this paper is the same as in |2^ and |2^, but different from one in 
]. This definition emphasizes the symmetry properties. Function u^{x;a,b,c;d) in (2.7)-(2.8) is 
obviously symmetric in b and c. Applying (III. 36) and Bailey's transform (III. 36) of [|12| once again 
we obtain 



u^{x; a, b, c; d) 

{abcq^+^.q^-^+'/d- g). 



{abcq',q^+'/d; g)^ 

X ^abcq''~^;q~'',abcdq''^^,aq'',bq'',cq''; g, 

{abq", acq", bcq", abcdq''-\ q^-^+^/d, q^-"-'/d; q)^ 
{ab, ac, be, abcdq^^'^, q^~^^/d, q^~^/d; g)^ 



(2.10) 



(2.11) 
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/ q-'',q^~''/ad,q^-^/bd,q^-^/cd 

X m ; q, q 

\ q'^'^^/abcd, q^-^+'/d, q^-^-^/d 

_^ (g-^ q^-'/ad, q^-" /bd, q^-'/cd, abcq''+', abcq"-'; q)^ 
{ab, ac, be, q^~'^''/ abed, q^~^^/ d, q^~^/ d; q)^ 

(abedq"~^ , abq^, aeq" , beq^ 
; 9, q 
abedq^^ , abeq'^'^^, abeq^~^ 

These representations show that Uy{x; a, b, c; d) is actually symmetric in a, b, and c. Here we have 
used the standard notation, 



2r+2W^2r+l Ctl, 02, • • • , a2r-i; Q', t) 



2r+2^2r+l 



^ a, q^/a, -qy/a, ai, 02, • • • , a2r-i 

; q, t 

\ y/a, -y/a, qa/ ai, qa/ 02, . . . , qa/ a2r-i 



(2.12) 



for a very-well-poised basic hypergeometric series defined before. 

Let us discuss analyticity properties of the function u^{z) defined in (2.7)-(2.8). One can easily 
see, that for integers u = n = 0,1,2,... this function is just a multiple of the Askey-Wilson 
polynomial (1.1) of the n-th degree, 

(•_l)"g-n("-l)/2 

Unix; a, b, c; d) = —— r — rf"" p„(x; a, b, e, d), (2.13) 

[ab, ae, be; q)n 

where x = (g^ + q~^) /2. In this case the Askey-Wilson polynomials p„(a;; a, b, e, d) are symmetric 
with respect to a permutation of all four parameters a, b, e, and d due to Sears' transformation P|. 

On the other hand, if v is not an integer, the essential poles of the gy^y-solution in (2.7) coincide 
with the simple poles of the infinite product 

U+^/d,q'-^/d;qy\ 

Therefore, the function 

Vy{z) = Vu{x{z);a,b,e;d) := iq^'^^ /d, q^~^ /d; q) u„{x{z);a,b,e;d), (2-14) 

V / 00 

where Uu{x; a, b, e; d) is defined in (2.7)-(2.8), is an entire function in the complex 2;-plane. 

Let us mention also that function u,y{x;a,b,e;d) satisfies the simple difference-differentiation 
formula 

S . . ^ , -^ 2g (1 - q-n (1 - abedq^-') 

Uu[x[z);a,b,e;d) =— — - — — — -— — (2.15) 



6x{z) ^ ^' ' ' ' ^ {I - q)d {I - ab){l - ae){l - be) 

x^z,_i (x{z);aq^'\bq^'\eq^/^;dq^'^) 

where 5f{z) = f{z + 1/2) - f{z - 1/2) and x{z) = (g^ + q~') /2. 
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3. Solutions of Pearson Equation 



In order to rewrite equation (2.1) for the function (2.7)-(2.8) in self-adjoint form (2.5), we have to 

find a solution of the Pearson-type equation (2.6). In the case of the g-quaratic grid ^ — 2 (^^ ~'~ ^ 
this equation can be rewritten in the form 

p{z + l) ai-z) 



p{z) (r{z + l] 



(3.1) 



(1 - ag-^-i)(l - bq-'-^){l - cq-'-^){l - q'+^/d) 
It is easy to check that 

fo, ,„(,)^ fa".^-"^f~ , (3.2) 

Po{z) q'-q-' 

Pa{z + 1) 



Pa{z + 1) 1 - 



Pa{z) 1-aq 



q ^ for p„(z) = (ag^ag ^g^+7a,g^ 7a;g)oo, (3.3) 
for p„(z) = (ag^,ag"^;g)oo. (3.4) 



(See, for example, pS), p9[, and ||30| for methods of solving the Pearson equation.) Therefore, one 
can choose the following solution of (3.1), 

^ (g- - q-T\q'\ q'-^^d, q^-^d; q)^ 

{aq"", aq-"", q^+'^/a, q^~^/a, aq^, aq-^, bq'^, bq-"", cq"", cq"^] g)oo' 



where a is an arbitrary additional parameter. It was shown in ||3^ that this solution satisfies the 
correct boundary conditions for the second order divided-difference Askey-Wilson operator (2.5) 
for certain values of the parameter a. 

Special cases a = (i or a = g/ci of (3.5) result in the weight function for the Askey-Wilson 
polynomials (cf. 0, 0). 

One can rewrite (3.1) as 

p{z + 1) (1 - ag^)(l - 6g^)(l - g"Vc)(l - q-'/d) 



p{z) (1 - ag-^-i)(l - 6g— i)(l - q'+^/c){l - q'+^/d)' 
We shall use the corresponding solution. 



(3.6) 



^ (g^+Vc,g^-Vc,g^+-Mg^-7d;g)^ ^ ^^^^^ 
{aq',aq-',bq',bq-^;q)^ 



in Section 7. 



4. Orthogonality Property 



Now we can prove the orthogonality relation of the gy^y-f unctions (2.7) with respect to the weight 
function (3.5) established in |Q. Let us apply the following g- version of the Sturm-Liouville 
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procedure (cf. 0, |T0|, and [^]). Consider the difference equations for the functions Uij{z) 
u^{x{z); a, b, c; d) and u^{z) = Ufj_{x{z); a, b, c; d) in self-adjoint form, 



A 



A 



a{z)p{z^ 



a{z)p{z) 



Vui,{zy 



+ Xf,p{z) u^{z) = 0, 



+ Xup{z) Uu{z) = 0, 



(4.1) 



(4.2) 



Vxi(z) y \ /f-^ / \/x{z) 

where the eigenvalues A = A;^ and A' = A^ are defined by (2.4). Let us multiply the first equation 
by Ui,{z), the second one by u^{z), and subtract the second equality from the first one. As a result 
we get 

(A;, - A^) u^{z) Uy{z) p{z) Vxi{z) 
= ^[a{z) p{z)W {u^{z),u,{z))] 



(4.3) 



where 



W{u^{z), Uu{z)) 



U^{Z) Uy{z) 

VUf,{z) Vuu{z) 



Vx(z) Vx(z 
Vuiy(z) 



(4.4) 



""^^^^ Vx{z) -^^'-^ Vx{z) 
Uy{z) u^{z - 1) - u^{z) u^{z - 1) 
x{z) — x{z — 1) 



is the analog of the Wronskian 25 



We need to know the pole structure of the analog of the Wronskian W{u^,Ui,) in (4.4). Let us 
transform the m's to the entire functions f 's by (2.14), 

Ueiz) = Lp{z) Vs{z), (4.5) 

where £ = /x, z/ and 

^{z) = (q'+yd,q'-yd;qy\ (4.6) 

Thus, 

W{u^{z),u,{z)) = ifiz) ip{z - 1) W{v,{z),v,{z)), (4.7) 

where the new 'Wronskian", W{v^{z),v^{z)), is clearly an entire function in z. 

Integrating (4.3) over the contour C indicated in the Figure; where the variable z is such that 
z = i6/\ogq~^ and —it < 6 < ti (we shall assume that < g < 1 throughout this work); gives 

(A;, - A^) J^u^{z)ui,{z) p{z)Vxi{z) dz (4.8) 

= [ A [a{z) p{z) ^{z) ifiz -1)W {v^{z),v,{z))] dz. 

All poles of the integrand in the right side of (4.8) coincide with the simple poles of the function 

a{z)p{z)^{z)ip{z-l) (4.9) 

^ C;(g2z^gl-22.^^^ 

(ag^ aq-', q^+'/a, q^-'/a; q)^ 



10 



SERGEI K. SUSLOV 
Im z 




Figure 

X fag^, a(^~^ , bq'', bq^"'', cq^, cq^~^, q^^^^d, q^~^ Id] q) ^ . 

V / oo 

The integrand in the right side of (4.8) has the natural purely imaginary period T — 2m /log q when 
< g < 1, so this integral is equal to 

[a{z) i^{z) ^{z - 1) p{z) W {v^{z)Mz))] dz, (4.10) 

where D is the boundary of the rectangle on the Figure oriented counterclockwise. 

The analog of the Wronskian W{v^, v^) is an entire function. Thus, when max (|a|, |c|, \q/d\) < 
1, the poles of the integrand in (4.10) inside the rectangle in the Figure are just the simple poles of 
p{z) ai z = aQ and z = \ — q.q, where g""^ = a and < Re ao < V^- By Cauchy's theorem, 

^J^[<j{z)^{z)^{z-l)p{z)W{v,{z)M^))] dz (4.11) 



where 
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= Res/(z)U^ + Res 

f{z) =a{z) p{z) if{z) ^{z -l)W {v,{z)M^)) (4-12) 

(ag^ aq-^, q^+^ / a, q^-^ / a; q)^ 

W{v^{z)Mz)) 



X 



{aq^.aq^ ^.bq^.bq^ ^.cq^'jCq^ ^^q^+^/d^q"^ ^/d;q)^ 
Evaluation of the residues at these simple poles gives 



Res = Jnn^(^ - ao) f{z) 



logg (g, q, aa, qa/a, ab, qb/a, ac, qc/a, qa/d, q^/ad] q)^ 

and 

Res = _ lim_ - 1 + ao) f{z) 



(4.13) 



dW {v,{z)Mz))l=^. 



(4.14) 



log q {q, q, aa, qa/a, ab, qb/a, ac, qc/a, qa/d, q^/ad; q)^ 
But 

W{v,{z),v^{z))l^^^ = W{v,{z),v^{z))l^,_^^ (4.15) 

due to the last line in (4.4) and the symmetry Ve{z) — Ve{—z), x{z) — x{—z). Thus, the residues 
are equal and as a result we evaluate the integral in the right side of (4.8), 

a{z) ^{z) ^{z - 1) p{z) W {v^{z),v,{z))] dz (4.16) 

c 

AnidW (t^^(ao),^!/(«o)) 



log q (q, q, aa, qa/a, ab, qb/ a, ac, qc/ a, qa/d, q^/ ad; q)^ 
_ 47rzri {a/d,q/ad;q)^ W {u^{aQ),u^{aQ)) 
log q {q, q, aa, qa/a, ab, qb/ a, ac, qc/ a; q)^ 

by (4.7). 

Combining (4.8) and (4.16), we, finally, arrive at the main equation, 

(A;u - A^) / Ui^{z) u„{z) p{z) 5x{z) dz (4-17) 
J c 

Airid {a/d,q/ad;q)^ 
log q {q, q, aa, qa/a, ab, qb/a, ac, qc/a; q)^^, 
X W {u^{ao),Ui,{ao)) , 

where max(|a|, |6|, |c|, \q/d\) < 1 and < Re ckq < 1/2, = a. 

It is worth mentioning the important special case first. If both of the "degree" parameters p and 
u are nonnegative integers: p — m — 0,1,2, . . . and u — n — 0,1,2, . . . ; equations (4.17) and (2.13) 
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result in the following real integral 

Pm{cos 9; a, b, c; d) Pn(cos 6; a, b, c; d) 



L 



(4.18) 



de 



L2ie^ e-^iS^qeiO/d,qe-'yd;q) 
(ae*^, ae~*^, qe^^/a, qe~^^/a; q)^ 

^ {a/d,q/ad; q)^ 

{q, q, aa, qa/a, ab, qb/a, ac, qc/a] q)^ 

^ -A'Kq^/'^d W ipm {m a, b, c; d) , Pn {m a, b, c; d)) 

involving the Askey-Wilson polynomials. Here we use the notation 

r) = x{ao) = ^(a + a'^^ . (4.19) 

One can easily see that when a = d or a = q/d our equation (4.18) implies the orthogonality 
relation for the Askey-Wilson polynomials (1.7). 

The new orthogonality property for the g^y-function appears when both of the parameters and 
u are not nonnegative integers. It is convenient in this case to rewrite (4.17) in terms of the entire 
functions t^'s as follows 



/ 

Jo 



i'^(cos 9; a, b, c; d) v,^{cos d; a, b, c; d) 



(ae^^ ae-^^ be'^, be^'^, ce^^ ce^'^, qe'^d, qe-'^d; q)^ 

(ae*^, ae"'^, qe^^ /a, ge~*^/a; q)^ 

= (q, q, aa, qa/a, ab, qb/ a, ac, qc/ a, qa/d, q^/ ad; q) 

—Anq^/'^d W {v^ {rj; a, b, c; d) , {rj; a, b, c; d)) 
1 - g - A,, 

The hmiting case // — > of (4.20) is also of interest. From (4.4), 

- A^ 



(4.20) 



1^^^ A^ - \y 

dv^{z) Vv^{z) d {Vv^{z) 

V,AZ] 



^. dfi Wxi^z) d/j, \ Vx{z) 

dvy{z) Vvy{z) d (Vv^{z)\ 

-Vv{z) 



dv Vx[z) dv \ Vx{z) 
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Therefore, 



/■TT 

Jo 



{vi,{cos6; a, b, c; d)y 



de 



(4.21) 



(ae*^, ae"*^, ge*^/ a, qe~^^ / a; g)^ 
fg, g, aa, qa/a, ab, qb/a, ac, qc/a, qa/d, q^/ad; q) 

V / oo 



X 



Anq^/^d 

1-g 



Vi,{r);a,b,c;d,) [ — Vi,{x;a,b,c; d) 

V 37 



x=r) 



- i;^(77; a, 6, c; d) — v^(x; a, 6, c; d) 



x=r]. 



Finally, choosing the parameters \x and v as ^-solutions of the equation 

Q {a + a"^) ; a, 6, c; = 0, 



(4.22) 



we arrive from (4.20) and (4.21) at the orthogonality relation of the sV^y-functions under considera- 
tion. 



rn 

Jo 



v^(cos 9; a, b, c; d) v^{cos 9; a, b, c; d) 
{ae^^, ae-^^, be'^, be'^^, ce'^, ce-'^, qe^^d, qe-^^d; q)^ 



(4.23) 



X 



-2ie. 



[ae^^, ae"'^, ge*^/ a, qe~^^/ a; q) 



d9^0 



if H ^ u, and 



/"TT 

Jo 



{vu{cos9; a, b, c; d)y 



ae^^ ae-*^ 6e*^ be-^^, ce^^ ce-*^ ge^Vc^, qe-^yd; q)^ 

X : : ^ (19 

(ae*^, ae"^^, qe'^^/a, qe~'^^/a; q)^ 
— (q, g, aa, qa/a, ab, qb/a, ac, qc/a, qa/d, (f'/ad; g) 

( Vi,{x; a, b, c; d) 



(4.24) 



X — i -wr- Mv; a, b, c; d) 

1 - g dX„ 



x=ri 



if 11 = u, respectively. We remind the reader that max(|a|, |&|, |c|, \q/d\) < 1. This will be assumed 
throughout this work. 

One can easily see that Vs{x; a, b, c; d) with e = fi,^ in (4.23)-(4.24) are real-valued functions of 
X for real e which are orthogonal with respect to a positive weight function when all parameters a, 
b, c, d, and a are real, or when any two of a, b, c are complex conjugate and all other parameters 
are real. We shall usually assume that as well. 
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5. Properties of Zeros and Asymptotics 



In the previous section we have proved that the 8</?7-functions (2.7)-(2.8) are orthogonal if the 
"boundary" condition (4.22) is satisfied. Let us discuss now some properties of z/-zeros of the 
corresponding "boundary" function, 



v (« + « ^) ; a, c; d 



X 8</27 



(gi ''a/d,bc;q)^ 
91 



/ aq-'' 



I aq 



d 



d 



bd cd 



aa, a/a 



q, hcq^ 



I aq 



\ 



V ' V 

{qa I (i, q^~^ I ad^ ahcq^ / a; q) 
{abc/a;q)^ 



aq aq q^ "a q^ ^ 

——, — , ab, ac, ——, — - 
d d d ad 



( abc I abc 



X 8^7 



qa ' V q^^'^a ' 



'',abcdq'' 
g^/^a a a a 



I abc abc q^abc q 



l-u 



a 



V 



qa 



qa a ad 



ab, ac, be 



J 



(5.1) 



(5.2) 



We have used (III. 23) of |[T^ to transform (5.1) to (5.2), compare also (2.10). Again, v^{ri; a, b, c; d) 
is real-valued function of u when all parameters a, b, c, d, and a are real, or when any two of a, b, 
c are complex conjugate and all other parameters are real. 

Main properties of zeros of the function (5.1)-(5.2) can be investigated by using the same methods 
as in |]TT|, [|T^, |TB| and [|T^]. The first property is that under certain conditions the real- valued 
function Vy{ri;a,b,c;d) has an infinity of positive z/-zeros. In order to establish this fact, one can 
consider the large z/-asymptotics of the 8V^7-function in (5.2), 



/ abc 



8^7 



abc 



qa ' V g^/^a ' 



abc _^ ^ u-i <^ c 
q ,abcdq 

a a a 



\ 



q^l'^a 



I abc abc q^abc q 



qa 



qa a ad 



■, ab, ac, be 



I abc 



6</35 



qa V q ' (y- 



abc 



abc abc 



1 /9 ' ' ' 

q^i^a a a a 



I abc 



abc 

\ qa' \ qa 
aa, ab, ac, abc/ a; q) 
{ab, ac, be, a^; q)^ 



, ab, ac, be 



a 



q, a 



(5.3) 
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by (11.20) of [0. Therefore, 

. (1 J-;^;.;^.!^^;;- (5.4, 

X (gi-7arf;g)^ [l + o(l)], 
as z/ — ^ oo. But for the positive values of q and ad the function 

(q"-" /ad] qj^ 



oscillates and has an infinity of real zeros as u approaches infinity (see |T^ and ||TB|). Indeed, 
consider the points u = Un = ujq + n, such that q'^° = /?, where n = 0, 1, 2, . . . and q < P < 1, as 
test points. Then, by using (1.8) of |1T2| , 



a. t.c;d)= '7f-°f-^°-;'';^'^ (5.5) 



' oo 



x(-i)",""('^^i)/2l^^^ [1 + 0(1)], 

as n ^ oo, and one can see that the right side of (5.5) changes sign infinitely many times at the 
test points v = lOn as, v approaches infinity. 

Thus we have established the following theorem. 

Theorem 5.1. The real-valued junction Vy{ri;a,b,c;d) defined by (5.2) has an infinity of positive 
V -zeros when aq < d, a^ < 1, ad > {), and < g < 1. Also, all parameters a, b, c are real, or when 
any of two of a, b, c are complex conjugate and the remainder parameter is real. 

Now we can prove the next result. 

Theorem 5.2. Function Vij{ri:ia,b,c\d) defined in (5.2) has only real v-zeros under the following 
conditions: 

1. parameters a and d are real; parameters a, b, c are all real, or, if any two of them are complex 
conjugate, the third parameter is real; 

2. the following inequalities holds: 

(a) max{\a\, |c|, \q/d\) < 1, q^^"^ < a < 1; 

(b) ad > q, aabc < q. 

Proof. Suppose that uq is a zero of function (5.2) which is not real. Let ui be the complex number 
conjugate to uq, so that ui is also a zero of (5.2) because this function is real under the hypotheses 
of the theorem. 

Consider equation (4.20) with /i = z/q and u = Ui. The integral on the left does not equal zero 
due to the positivity of the integrand, but the analog of the Wronskian on the right is zero in view 
of (4.22). Therefore, 

= ^fi, (5-6) 
the eigenvalues defined by (2.4) are real. The last equation can be rewritten as 

(l - g'^^-'^") (l - abcdq'"'+'''-^) = 0. (5.7) 
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The first solution is uq = ui, so the roots are real in this case. 
The second solution of (5.7) is 

-.1^1 



q 



q 



q 



abed ' V ahcd 

where ahc > and x is an arbitrary real number. But in this case our function can be represented 
as a multiple of a positive function. Indeed, when u = uq or u = ui, 

v^{r]]a,b, c; d) 

{q/ad, q^~"a/d, aahcq^; q)^ 



(5.9) 



{aahc; q)^ 

I aabc laabc I aabc , , 

y "^172"' ~y "^172"'^ ,a6cdg ,aa,ab,ac 



X s^7 



\ 



/ aabc / aabc 



q_ 

ad 



\ 



q 



q 



q'^aabc, q "a/d, ab, ac, be 



J 



{q/ad,q ^a/d,aabcq^]q)^ 



X 



1 + E 



[aabc] q)^ 

1 — aabcq^'^~^) {aabcq~^, aa, ab, ac; q)^ 



n=l 



{1 — aabcq ^) {q,ab,ac,bc]q)^ 
{q^" ,abcdq^^^]q)^ ( q \^ 



X 



{q'^aabc, q^~'^a/d; q)^ \ad 



by (5.2), (III. 23) and (1.31) of |T^. One can easily see that under the hypotheses of the theorem 
all products in the last but one line of (5.9) are positive. From (1.3)-(1.4) and (5.8), 

n-l / 



I abed u abed ,1, 

cos xq + 



k=0 
n-l 



q 



q 



q 



fc=0 



I abed 



abed 



2 „2fc 



and 



n— 1 



{^q^aabe,q^ "a/d^q^ = Y\ j 1 — 2a;^ 



cos xq \ ^ sm X q 

y q 



labeq ,^ ^abeq 

cosx q + tt — T' q 



(5.10) 



k=0 
n-l 

n 

k=0 



d 



d 



I — a\ 



labeq 



V 



d 



COS X q + «■ 



, abeq . 



d 



2 ^2k 



sm xq 



(5.11) 



where n = 1, 2, 3, . . . , 00. Thus, these products are positive too, which proves the positivity of the 
"boundary" function in (5.9). So we have obtained a contradiction and the complex zeros (5.8) 
cannot exist. This completes the proof of the theorem. □ 

Theorem 5.3. The real v-zeros of v^{ri]a,b,e;d) are simple under the hypotheses 1 and 2(a) of 
Theorem 5.2. 
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Proof. This follows directly from the relation (4.21). Indeed, the integral on the left side is positive 

d 

under the conditions of the theorem, which means that — fy(?7; a, b, c;d)j^O when v^{ri; a, b, c; d) = 

0. ^ n 

Let us consider two functions, 

f\u) = Vy{ri]a,b,c]d), (5.12) 

(5.13) 



V 

g{i') = —v^{x;a,b,c]d) 
vx 



x=rj 

Our next property is that positive zeros of /(z/) and (?(//) are interlaced. 

Theorem 5.4. // Ui, z/2, 1^3, ■■ ■ are the positive zeros of f{v) arranged in ascending order of mag- 
nitude, and /ii,/i2,/i3, . . . are those of giy), then 

< /ii < Z/i < /i2 < z/2 < /is < ^^3 < • • • , (5.14) 
if the conditions of Theorem 5. 3 are satisfied. 

d 

Proof. Suppose that Vk and vt+i are two successive zeros of /(z/). Then the derivative 'g~fi^) has 

different signs at z/ = z/^ and u = Vk+i- This means, in view of the positivity of the integral on the 
left side of (4.21), that g^v) changes its sign between z/^, and Vk+i and, therefore, has at least one 
zero on each interval (z/^, Vk+i)- 

To complete the proof of the theorem, we have to show that g{v) changes its sign on each interval 
{i^k,i^k+i) only once. Suppose that g{iik) = g{t^k+i) = and Uk < fik < f^k+i < ^k+i- Then, by 
(4.21), function /(z/) has different signs at z/ = /ifc and z/ = fik+i and, therefore, this function has 
at least one more zero on (z/^, z/^+i). So, we have obtained a contradiction, and, therefore, function 
g{i') has exactly one zero between any two successive zeros of /(z/). 

In order to finish the proof of our theorem one has to show that fii < ui. We would like to leave 
the proof of this fact as a congecture for the reader. □ 

Let us discuss an asymptotic behavior of u^{x] a, b, c; d) for large values of z/ applying the methods 
used in |^T|, and [jT^ at the level of the Askey- Wilson polynomials. In a similar manner, the 
8V?7-function in (2.10) can be transformed by (III. 37) of ||T^, 

sWr (abcq"^; abcdq"-^ , aq\ bq\ cq'; q, ^ j (5.15) 



^,bq ^,cq ^ ,bq^^^^^ , cq^^^^^ , dq^ ^ ,abcq^ ^bcdq^^^-^q)^ 



+ 



(g~2^, afe, ac, 6c, q^^^ , bdq^, cdq'^, 6cg^+^+^^; q)^ 

X sWj {bcq"^^'] bcq", bq\ cq\ q^+'/a, q^+'/d; q, adq") 

{aq^, bq^, cq^, bq^^^~^, cq^^"'^ , dq^^^, abcq^, bcdq"'^] q)^ 
(g^^, ab, ac, be, q^^" , bdq'^, cdq^ , bcq^^"^"^""] q)^ 
{q^-'-'/d,q^+'/d,abcq^-'-q)^ 
[q^-'^/d, q^-^+^/d, abcq''+''; q)^ 

X sWj{bcq''-'^']bcq''M'\cq-',q^-'/a,q^-'/d]q, adq") 
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Therefore, 

u^{x]a,b,c;d) (5.16) 
_ (6gi+'^+^ cgl+^+^ dq"-', ahcq"^^, bcdq''+'; q)^ 
~ {q^+'',bdq'',cdq'',bcq^+''+^'^;q)^ 
^ {aq'\ bq~^, cq~\ q^-^^'/d; q)^ 
{ah, ac, bc,q-^',q^+'/d; q)^ 

X gW^T {hcq"^^'- bcq", hq\ cq\ q^^' j a, q^^' j d' g, adq") 

{bq^+^-', cgl+'^-^ abcq''-', bcdq"-'; q)^ 

{q'^+'',bdq'',cdq'',bcq^+''-^'';q)^ 
^ {aq',bq',cq',q^-''-'/d] q)^ 
{ab, ac, be, q^'^ , q'^-' /d; q)^ 

X ^Wj {bcq"-^'; bcq", bq-\ cq'', q^-'/a, q^-'/d; q, adq"") . 

We shall use the last expression to determine the large u asymptotic of our orthogonal s 927- function 
later. 

For the "boundary" function (5.2) equation (5.16) takes the form 

■i^i/ Q (« + a""^) ; a, c; c?^ (5.17) 

{abq^^''', acq^^" , d(f j a, aabccf , abcdq'^; q)^ 
~ {q^+'',bdq'',cdq'',a'^bcq'^+'';q)^ 
{a/a, b/a, c/a, q/ad, aq^^" /d; q)^ 
{ab, ac, be, a'"^; q)^ 

X ^Wj {a%cq'^]bcq'' ,ab,ac,qa/ a,qa/ d]q, adq" 
[Pq^+^ I Q,^ cq^+'' /a, adq^, abcq^ / a, bcdq" / a; q)^ 
{q^+^,bdq^,cdq'',bcq^+'' /a^]q)^ 
{aa, ab, ac, qa/d, q^~^ jad; q)^ 



X 



X 



{ab, ac, be, a"^; q)^ 

X gWj (beq'^ / o?\ beq'^ , b/ a, cj a, qj aa, qj ad; q, adq'^^ . 

The last term dominates here as u approaches infinity (cf. (5.4)). 

The proof of Theorem 5.1 has strongly indicated that asymptotically the large positive z/-zeros 
of Vu{ri; a, b, c; d) are 

iyn = n + en, < e„ < 1 (5.18) 

as n ^ 00. The same consideration as in [|n|, |T^, and |]1T| shows that this function changes sign 
only once between any two successive test points u = Un and u = cUn+i defined on the page 15 for 
sufficiently large values of n. We include details of this proof in Section 9 to make this work as 
self-contained as possible. 

Our next theorem provides a more accurate estimate for the distribution of the large positive 
zeros of the "boundary" function (5.2). 
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Theorem 5.5. // z/i, z/2, z^s, . . . are the positive zeros of Vy{rj] a, b, c; d) arranged in ascending order 
of magnitude, then 

\og(ad) , 

iyn = n- +0(1), 5.19 

as n —y 00. 

Proof. All zeros of our function (5.2) coincide with the zeros of a new function, 

w,^(r]] a,b, c; d) := (—q^^" /ad; q) Vi^(r]] a,b, c; d), (5.20) 

because {—q^~^/oid]q)^ is positive for real u. Equations (5.17) and (5.20) give us the large z/- 
asymptotic of w^irj] a, 6, c; d). When z/ = 7^ such that g"^" = q^ /ad and n = 1, 2, 3, ... , the second 
term in (5.17) vanishes and we get 

{qa^,a/a,bc,qb/d,qc/d,abc/d;q)^ 

w^n iv, a, b, c; d) = 5 21) 

(— g, ao, ac, oc, abcq/d; q)^ 



{q/ad, b/a, c/a, abcq/d; q)^ 
(—1, 6c, qc/d, abc/d; q)^ 



{---)' 



abc ^ b ^ a a a ^ 

X gVr7 — g";-g",a6,ac,g-,g-;g, - 

\ d a a d a , 



with the help of (1.9) of Thus, 

lim (n: a, b, c: d) = 0, (5.22) 
which proves our theorem. □ 

Let us discuss also the large positive //-asymptotic of Uy{x\a,b,c\d) when x = cos 9 belongs to 
the interval of orthogonality — 1 < x < 1. It is clear from (5.16) that the leading terms in the 
asymptotic expansion of this function are given by 

(ae'^,be''^,ce'^,q^-''e-'yd; q) 

w,(cos 6; a, 6, c; d) ~ ^— , /"^ (5.23) 

[ab, ac, be, e^*'', ge *'^/ d; q)^ 

(ae-'^, be-'\ cc'^\ q^'^e'^d; q) 
I _V / 00 



{ab, ac, be, e~^*^, qe'^^ /d; g)^ 

In particular, when u = Un are large zeros of Uult]; a, b, c; d), we can estimate 

'U;^,^(cos6'; a, b, c; d) ~ m^„(cos6'; a, b, c; (i), (5.24) 

where g^" = q"' /ad, due to (5.19) as n ^ 00. Relations (5.23)-(5.24) lead to the following theorem. 

Theorem 5.6. For —l<x = cos 6' < 1 and \q\ < 1 the leading term in the asymptotic expansion 
o/'U^,^ (cos 6*; a, b, c;d) as n ^ 00 is given by 

g-(«-l)("-2)/2 

(cos 9; a, b, c; d) ~ — — (5.25) 

{ab, ac, be; q)^ 



X 



A{f')\ cos((n-l)^-x), 
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where 



and 
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A(e 



(qe^^jd, qe-^^/d;q) 

V / CXD 

(ae*^, ae~^^, qe'^^/a, qe~'^^/a] q)^ 



X = ar} 



(5.26) 
(5.27) 



(5.28) 



It is worth mentioning that (5.27) coinsides with the weight function in the orthogonahty relation 
for Ui, (cos 9;a,b,c;d). 

In a similar fashion, one can consider some other properties of zeros of the sV^r-function (5.1)- 
(5.2) close to those established in [|1^ and [|17| at the level of the basic trigonometric functions 
and the g-Bessel function, respectively. 



6. Evaluation of Some Constants 

In this section we shall find an explicit expression for the squared norm, 

^2 ^ r ul{cose;a,b, c; d) 

Jo (ae*^ ae-*^ be'^, be~'\ ce\ ce"*^; q)^ 



e-2*«, ge^Vrf, qe-'^jd; q 



^d9, 



(ae*^,ae ^^,qe^^ /a,qe ^^/a;q)^ 

on the right side of (4.21). It is more convinient to write the Mj,'s instead of v^^s here, cf. (2.14), 
(4.21), and (6.1). Using (2.8) with a and c interchanged we get 



2 {abq'',q^ '"/cd;q)^^ 



dt = ' ,' {'^ (6.2) 

(a6, q/cd] q)c 



loo 

^ - {q-'^,abcdq^-\qX 
^0 (g, ac, be, cd; g)„ 



X 



Uu{cos 9; a, b, c; d) 



y _V / CXD 

(ae*^, ae~*^, qe^^ /a, ge~*^/ a; g)^ 
_^ {q-" , abcdq''-\ qa/d, qb/d; q)^ 



{ab, ac, be, ed/q; q)oo 
(a6g^gl-7cd;g)„ 



n=0 



{q,qa/d,qb/d,q^/ed; q)n 
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u^{cos 9; a, b, c; d) 
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e ,qe /d, qe /d; q 



X 



Both integrals on the right side have the same structure and can be evaluated in a similar way. Let 

Mi,(cos6'; a, b, c; d) 



fe^^^ e-^'\qe'^/d,qe-'^/d] q) 
(ae*^, ae~*^, qe^^/a, qe~'^^/a; q)^ 
{qa/d,bcq'';q)^ 



(6.3) 



(gi ''a/d,bc;q)^ 



X 



a 



sVTt ^ q""; q^'Vbd, q^-^/cd, ae"*'; g, 6cg^ 



ae*^, ae~*^, 6e*^, fee~*^, 7g"'e*^, 7g"e~*^; g)^ 
(gi-^Vd, q^-''e-'yd;q)^ 



X 



(ae*^, ae~*^, qe^^/a, ge~*^/a; g)^ 



by (2.7) with 7 = c and = q/d. The last integral can be evaluated in terms of two balanced 5V24's 
by (6.4) of 13, 



£ sWr (dg/q- h, r, g/f, de'\ de-''; g, gf/hi 



(6.4) 



X 



(e2^^e-2^^;g 



-ie 



be-'^, ce-^\ de^ , de-'^\ g)^ 
(^e*^ ^e"*^; g] 



+ 



(/e^^ /e-^^; g)^ 
27r {abed, dg,g/d] q)oo 
(g, a6, ac, ad, be, bd, ed, df, f /d; g)oo 

^ ad,bd,ed, g/ f,dg/hr 
X 504 ; g, g 

abed,qd/ f,dg/h,dg/r 

2n {abef, dg,g/f;q)^ 
(g, ab, ae, af, be, bf, ef, df, d/ f; g)oo 
^ {.fg/h,fg/r, dg/hr;q)^ 
{dg/h,dg/r, fg/hr;q)^ 

( af,bf,ef,g/d,gf/hr 



X 504 



V abef,qf/d,gf/h,gf/r 



• q, q 
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this formula appears in a straightforward manner if one uses Bailey's transform (III.36) of and 
the Askey-Wilson integral. Thus, 

27r [qa/ d, q^'" / ad, bcq", a'jq; q)^ 



(g, q, ab, be, b/ a, aa, qa/ a, 07, q'j/a, 07; q)^ 
(a7,g7/a,a7;g)„ 



(6.5) 



X 



X 504 



+ 



f a'yq"',acq",q^~" /bd,aa,qa/a 

\ a'yq^^'^,ac,qa/b,qa/d 
2tt {qb/d, q^-^/bd, acq", kyq; q)^ 
(g, q, ab, ac, a/b, ab, qbj a, a^, q^ /cy., 67; q)^ 
{a-f,q-f/a,b-f;q)^ 
iPiq; q)n 

( hyq^,bcq" ,q^~" /ad,ab,qb/a 



X 504 



g, Q 



\ 679^+", be, qb/a,qb/d 

where 7 = c and 7 = q/d. One can see that the second term here equals the first one with a and b 
interchanged. 

Combining (6.2), (6.3), and (6.5), we obtain 

f^(cos6'; a, b, c; d) 



d. 



(ae^^ ae-^^ be^^ , be'^^ , ce^ , ce-'^, qe'^d, qe-^yd; q)^ 



(6.6) 



■dd 



27r {qa/d, abq" , bcq" , q^~"/ad, q^~"/cd; q)^ 



(1 — ac){q, ab; qY (be, b/a, q/cd, aa, qa/a, ac, qc/a; q)^ 
„ {q-", abcdq"^^, ac, qcj a; g)„ 

n=0 



00 



(g, acq, be, cd; q)n 
( acq", acq", q^~"/bd,aa, qa/a 



X 504 



+ 



y acq"^^ ,ac,qa/b,qa/ d 
277 [qb/d, abq", acq", q^-" /bd, q^'" jcd; q)^ 



1 — bc){q, ab; g)^ (ac, a/b, q/cd, ab, qb/a, ac, qc/a; g)^ 
{q~", abcdq"~^, ac, qc/ a; g)„ 

J ■ 

n=0 



00 

X 



(g, qac, ac, cd; g)„ 
/ bcq^,bcq" ,q^~" /ad,ab,qb/a 



X 504 



+ 



y bcq"''^^, be, qb/a, qb/d 
2tt {qa/d; q)1^ {q''^, abedq"^^, q^^"/ad, bcq", qb/d; g)^ 



^1 — qa/d){q, ab, be; g)^ (ac, b/a, cd/q, aa, qa/a, qa/d, q^ /ad; q)^ 
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^ y n Vcd, qa/d, q'^/ad; q)^ 

^0 {q, Qb/d, q^a/d, q^/cd; q)^ 

^ aq^^^ / d, acq'' , q^"'' /hd, aa, qa/a \ 
X 504 ; q, q 

y aq^^"^ /d,ac,qa/h,qald j 

^ 2-71 {qb / d; q)1^ {q'" jabcdq"^^ jq^^" /bd, acq" , qa/d; q)^ 



1 — qb/d){q, ab, ac; qY {be, a/b, cd/q, ab, qb/a, qa/d, q^ /ad; q)^ 
{abq'^, q^^" / cd, qa/d, q^ / ad; q)^ 
n=o {q,qa/d,q^b/d,q^/cd;q)n 



oo 

X 



bq"'^^ /d,bcq'^ ,q^ / ad,ab,qb/ a 

X 504 I ; g, g 

bq^'+^/d, be, qb/a, qb/d 

One can see again that the second and the forth terms in this formula are equal to the first and 
the third ones, respectively, with a and b interchanged. When u satisfies the boundary condition 

(4.22) the last integral gives the values of the normalization constants in the orthogonality relation 

(4.23) -(4.24). 

7. Some Identity 

In this section we shall derive an interesting relation involving a determinant of four gf/^y-functions 
of type (2.7). Let 

u{z) = Uu{x{z);a,b,c;d), (7.1) 

v{z) = u^{x{z); a, b, d; c) = u{z)\^^^ (7.2) 

be two solutions of equation (2.1) corresponding to the same eigenvalue (2.4). Then, due to (4.3), 

A[a{z)p{z)W{u{z), ^;(^))]=0, (7.3) 

where W{u,v) is the analog of the Wronskian defined in (4.4) and p{z) is the appropriate solution 
of the Pearson equation (3.7). Using the difference-differentiation formula (2.15) we can rewrite the 
"Wronskian" as 

T.. / / N / ^^ 2g (1 - q-") (1 - abcdq"-') 

W iuiz), viz))^- \t7^ , (7.4) 

^ ^ ^" {l-q)c{l-ab){l-ad){l-bd) ^ ' 

X (x{z - 1/2); agV2, hq^l\ dq^'^; cq^''') 
X Ui, {x{z); a, b, c; d) 

2q (1 - q-") (1 - abcdq"-^) 
~ {l-q)d {1 - ab){l - ac){l - be) 

X u^_, {x{z - 1/2); agV2, bq^^^ cq'l^; dq"') 

X Uv {x{z); a, b, d; c) . 

One can easily see that the function 

g{z)^a{z)p{z)W{u{z), v{z)) (7.5) 
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(g7c, g^'Vc, gVd, q^-^/d\ q)^ 



{aq\aq^-\hq\hq^-^]q)^ 
X W{u{z), v{z)) 

in (7.3) is doubly periodic function in z without poles in the rectangle on the Figure. Therefore, 
this function is just a constant by Liouville's theorem, 

a{z) p{z) W {u{z), v{z)) = C. (7.6) 

To find the value of this constant we can choose here z = zq such that = a. From (7.4) and 
(2.7) one gets 

W{u{z^), v{zo)) 
_ 2g(l - g-^) (1 - abcdq^-^) 

(g, qbc, qb/d, qc/ d, bdq^ , (^~^ I ac, q^~^a/d, a?bcq"] q)^ 



(7.7) 



X 



acd {ab, qac, be, bd, 1/ac, q/ad, a/d, qabc/d; q)^ 

X (^abc/d;ab,ac,a/d,q^~'^/ad,bcq'^;q, q^ 

(g, qbd, qb/ c, qd/c, bcq'^, q^~^ / ad, q^~^a/ c, o^bdq^; q)^ 
acd {ab, qad, bd, be, 1/ad, q/ac, a/c, qabd/c; q)^ 

X (abd/c; ab,ad,a/c,q^~'^ /ac,bdq'^; q, qj 

2q (g, a^, qb/a, d/c, qc/d, q~^ , abcdq^^^, abq^ , q^^^ /cd; g)^ 
(1 — q)acd {ab, qac, ad, be, bd, 1/ac, q/ad, a/c, a/d; g)^ 
by (III.24) and (11.25) of [0. 

As a result, from (7.6) and (7.7) we find the value of the "Wronskian" of the sV^y-f unctions (7.1) 
and (7.2), 

2g {c/d, qd/c, q~^ , abcdq^~^ , abq^ , q^~^ /cd; g)^ 



W{u{z), v{z)) 



X 



(1 — g)c [ab, ab, ac, ad, be, bd; q)^ 
{aq\aq^~\bq\bq^-'';q)^ 



(7.8) 



{q^/c,q^ ^/c,q^/d,q'' '/d;q)^' 
Due to (7.4) and (2.7) the last equation can also be rewritten in a more explicit form, 

d {ac, adq; g)^ 
(q^-a/cq^-'^a/d; g)^ 

(g^-^+Vc, g'-'^-Vc, q'-'+^'/d, q'-'^-^d; g)^ 



(7.9) 



(gVc,g^-Vc, q'+'/d,q^-'/d; g)^ 
X {q'"a/c; q''\ q^'^ /be, q'-'/cd, aq\ aq'-'; q, g) 

X (q^" a/d; q~'',q^~''/bd,q^~''/cd,aq\aq~'';q, g) 

_ c{ad,acq; g)^ 
{q^-''a/d,q^-''a/c; g)^ 



X 
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(gi-'^+Vrf, g^-'-yd, g^-^+Vc g^^-yc; g)^ 
{q-/d,q'--/d,q'+^/c,q'--/c- g)^ 

sW^T {q^'^a/d; q^-\ q^'^/bd, q^-^/cd, aq\ aq^-'; q 
X (^q~'^a/c;q^'^,q^^'^ /bc,q^^'^ /cd,aq'',aq^'';q, q 

(c/ d, qd/c, q^~^ , abcdq" , ahq" , q^~^ / cd; g)^ 

c (a6, a6g, qa/c, qa/d, bcq", bdq"; q)^ 

^ d{aq'',aq^~'',bq'',bq'^~'';q)^ 

{q'/c, q^-'/c, gy d, g^-y rf; g)^ ' 

The second term on the left side here is the same as the first one with c and d interchanged. 

One can see from (7.8) that two solutions u{z) and v{z) are linear dependent when u is an integer. 
In this case due to (2.13) both solutions are the Askey-Wilson polynomials, up to a factor, which 
are related by Sears' transformation. On the other hand, equation (7.8) shows that there is no 
analog of Sears's transformation at the level of very-well-poised 807-functions. 

8. Some Special and Limiting Cases 

The Askey-Wilson polynomials are known as the most general system of classical orthogonal 
polynomials. They include all other classical orthogonal polynomials as special and/or limiting 
cases ITII, [ 



fT^ , |23|, and [25]. Let us discuss in a similar manner a few interesting special cases 
of the orthogonal sV^r-functions (2.7). 

8.1. Extension of continuous dual g-Hahn polynomials. Letting c ^ in (2.7) and then 
changing by c we get 

(ga/c,gi-'^+yc,gi-'^-yc;g)^ 



u^{x; a, b; c) 



[q^-^a/ c, gi+y c, g^-y c; g)^ 



X jip7 



aq 



I aq 



I aq 



-u 1 



l-u 



\ 



be 



,aq ,aq 



I aq 



V 



Uq-" aq , q^-""^^ q^-"'' 

, — , ab, 0, , 

c c c 



(g^ "/aqq)^ 
{q/ac] q)^ 



3V'2 



c 



g '',aq^,aq 



V 



ab, ac 



+ 



{q-",qb/c, aq'',aq-';q)^ 
{ab, ac/q,q^+^/c,q^-^/c; g), 

q^'"/ ac, g^+^/ c, g^~^/ c 

X 3V^2 ' 



qb/c, q^ /ac 

The 7937-function here can also be transformed to a 3^92 by (3.2.11) of |T^, 

(a6g^gl-'^+yc,gl-'^-yc;g)^ 



Uy{x; a, b; c) 



(a6,gi+yc,gi yc;g)^ 



(8.1) 



(8.2) 



26 SERGEI K. SUSLOV 

X 3V92 ; g, ahq"" \ . 

One can easily see that for an integer v our function (8.1)-(8.2) is just a multiple of the continuous 
dual g-Hahn polynomial [Q , , , and . 

The orthogonality relation for the corresponding entire function, 

Vu{x;a,b;c) = [q^^'' /c, q^'"" /c; q) Uu{x]a,b]c), 

V / oo 



takes the form 



^^^(cos 9] a, b] c) f ,^(cos 9; a, b; c) 



( 



8.3) 
8.4) 



if fi ^ u, and 



[ae*^, ae~*^, ge*^/ a, ge~*^/ a; g)^ 



(fi^(cos6'; a, b; c))^ 
(ae*^, ae~*^, 6e*^, 6e~*^, qe^^ /c, ge~*^/c; g)^ 



c/e = 



(8.5) 







(ae*^, ae~*^, qe^^/a, ge~*^/a; g) 
(^g, g, aa, qa/a, ab, qb/ a, qa/c, q^ / ac] q 
-A-nq^/^c d , , . / V 



■v„{r];a,b;c) (^—v„{x;a,b;c)^ 



x=ri 



1 - g dX, 

if fi = respectively. The "degree" parameters /i and u here satisfy the "boundary" condition 

(8.6) 

/ q^'^,a/a,b/a 



^(a + a ^) ; a, 6; cj 

(qa/c,q^-yac;q) 
V / 00 



a 

g, q 

A^e / I. C 



0. 



y g ^/ac,ab 

Properties of zeros of this function follow as a special case of the results of Section 5. 

8.2. Extension of Al-Salam and Chihara polynomials. Letting 6 ^ in (8.1) (8.2) and then 
changing c by 6 one gets 

(gi-'^+V6,gi-'^-V6;g)^ 



u^{x; a; b) 



iq^+^/b,q^--/b- g)^ 
/ q-^,q^-^/ab 



(8.7) 



\ q^+'/b,q^-'/b 
{q^-^/ab-q) 



g, g^ 



[q/ab] q)^ 



3V2 



V 



g ,aq,aq 
ab,0 



g, g 
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ah/q,q^+^/h,q^-''/h] g)^ 

q^-''/ab,q^+'/b,q^-'/b \ 

; q ■ 

g^/a6, y 

For an integer u this function is a multiple of the Al-Salam and Chihara polynomial IQ, 0], 0, and 

The orthogonality relation of the entire function, 

Vu{x;a]b) = (q^^yb,q^~yb;q) u^{x]a]b), (8.8) 



IS 



v^{cos9]a;b)vu{cos9]a]b) 
(ae*^, ae~*^, qe^^/b, qe~^^/b] q)^ 



(8.9) 



X 



= 0, 



where v are solutions of 



^ (a + a ^) ; a; 6 



(8.10) 



= [qa/b, q^ ^ /ab] q 
For properties of zeros see Section 5. 



/ g ^,a/a 



a 



8.3. Extension of continuous big g-Hermite polynomials. Letting a — in (8.7) and then 
changing 6 by a we have 

(gi-'^+Va,gi-'^-Va;g)^ 



Uy{x] a) 



{q^+^/a,q^ ^/a;q) 
( ? 



(8.11) 



oo 

V 



\ gi-'^+Va, gi-^-V 



For an integer z/ this function is a multiple of the continuous big g-Hermite polynomial |23]. 
The orthogonality relation has the form 

r-K Vfj_{cos6;a)v^{cos6;a) (^e^*^, e~^*^; g^ 



where ^ v satisfy 



(ae*^,ae ^^.qe^^ja^qe "^^ j a,qe^^ j a,qe ^^/a\q) 

1 / -1 
-\a + a 



^— de = 0, 



(8.12) 



(8.13) 



(g^ yaa;q) lif^ 



/ q/aa 
\ q^^^/aa 



g, g ^a/a = 
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and t;£(a;;a) = (g^+Va,g^ V'^;9)oo ^eix^a). 

8.4. Extension of continuous g-Hermite polynomials. The continuous g-Hermite polynomi- 
als Hn{x\q) are the simplest special case a = b = c = d = Oof the Askey- Wilson polynomials 
Pn{x;a,b,c,d) or the special case a = of the continuous big g-Hermite polynomials Pn{x;a) [H, 
^ ^ n. See H and § for the proof of 



T|], [ID, and 



linip„(x;a) = i?„(x|g) 
directly from the series representation of the continuous big g-Hermite polynomials. 

Let us consider the difference equation (2.1) with a = b = c = d = and let us choose the 
following solution, Uu{z) = {x{z)\q), such that 

HMq) 



-q 



X 2V52 



1+22 



(8.14) 



V 



q~ 



.\-v\2z 



q 



?^ q 



-q- -g^ 

as a nonterminating extension of the continuous g-Hermite polynomials. This solution differs from 
the corresponding one in by a periodic factor. For an integer v function (8.14) coinsides with 
Hn{x\q) up to a constant. 



Comparing (8.14) with the equation (2.8) of [rT| one can see that function Hy{x\q) is a multiple of 
the basic cosine function C{x] uj) for certain values of parameter v. Therefore, this function satisfies 
the orthogonality relation (1.16) under the "boundary" conditions (1.20). We would like to leave 
the details to the reader. 

In a similar fashion one can consider some other special and limiting cases of the orthogonal 
8 V97-f unctions. 



9. Appendix: Estimate of Number of Zeros 

In this section we give an estimate for number of zeros of the "boundary" function Vu{ri\ a, 6, c; d) 
on the basis of Jensen's theorem (see, for example, p^). We shall apply the method proposed by 
Mourad Ismail at the level of the third Jackson g-Bessel functions (see also and |11[] for 
an extension of his idea to g-Bessel functions on a g-quadratic grid and g-trigonometric functions, 
respectively). 



Let us consider the entire function 

/(C) = Q (a + a"^) ; a, 6, c; d 
[qa/d- q)^ 



(9.1) 



[ahcja; g) 



oo 



oo 



X n (1 - Cq'^^^/oid + ahcq 

k=0 

^ / a\"' n-abcq'^"'-^ 

X E (q^ 



m=0 



"^^'/a'd) 

/a) {abcq~^/a, a/a, b/a, c/a; q)^ 



1 — abcq^^/a) (g, ab, ac, be; g)^ 
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m—1 



X 



Hi 



1 — (q'^ + abcdq 



2k-l 



k=0 



in a complex variable 



C = q^" + abcdq"'^, 



(9.2) 



where \qa/d\ < 1 (cf. (5.2)). 

Let Uf (r) be the number of of zeros of f{Q in the circle \C,\ < r. Consider also circles of radius 
R = Rn = q~"'/P + [3ahcdq^~^, q < (3 < 1 with n = 1, 2, 3, . . . in the complex ^-plane. Since nj (r) 
is nondecreasing with r one can write 

nj (Rn) < rif (r) < nj (Rn+i) 

ii Rn ^ r < Rn+1, and, therefore, 



Rn r 



But 



R. 



n + l 



Rn 



Rn+i rif (r) 



logg ^ + o(l) 



, , q~^ + P'^abcdq^'' 
- = logr = log I 

'Rn r 

as n — i> oo, and, finally, one gets 

logg"^ Uf (Rn) < / 

J B 

for sufficiently large n. 

The next step is to estimate the integral in (9.6). By Jensen's theorem 

rRn+i Uf {r 



' Rn 



dr < log q ^ Uf [Rn+i] 



dr 



Rn r 



i^n+i rif (r) 



dr 



Rr 



Uf (r) 



dr 



r 



1 



27r Jo 



log 



/(i?ne 



For large values of n we have i?„ ~ and, in view of (5.3) and (9.1), 



n+i e 



fiRne" 



St 



1 - e^'^ q^-''/a(3d + ahcq^^^^ / d 



k=o - e'' q^-'^+^/al3d + abcq'^'^+^a^d 
-e'^q-^'/aPd, 



where we have used the formula 

oo 

n (l - + 9q 



k=0 



2k 



n-1 



fc=0 



k=n 



fc=0 



^ h h 



(9.3) 
(9.4) 

(9.5) 
(9.6) 



(9.7) 



(9. 



(9.9) 
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X 



n (l - hq' + 



2n+2fc 



fe=0 



with h — qe^'^ /a(5d and g — qabc/a^d. Thus, 



log 



where 7 = |Q;/3(i|, and 



^"+1 Uf (r) 



~ n log q ^ — log 7, 



Rn 



dr — nlogq ""^ — log7 + o (1) 



(9.10) 



as n ^ 00. 

Prom (9.6) and (9.10), 

^ _ log 7/ log - 1 < ^/ (^") < 1 _ log 7/ log 

and, therefore, 

hm = 1. 

n->-oo n 



n 



(9.11) 

(9.12) 
(9.13) 



On the other hand, from (9.11), 

n — 1 — log7/ logg~^ < n/ (i?„) < n — log7/logg~^. 

The difference between the upper and the lower bounds here is 1 which means that there is only 
one positive root of v^{ri;a,b,c;d) between the test points i/ = a;„ and u = uJn+i defined on the 
page 15 during the proof of Theorem 5.1 for large values of n. 
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